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Abstract 

We provide a new estimate of the long-distance component to the radiative 
transition B —* pj. Our mechanism involves the soft-scattering of on-shell 
hadronic products of nonleptonic B decay, as in the chain B — > pp — > pj. We 
employ a phenomenological fit to scattering data to estimate the effect. The 
specific intermediate states considered here modify the B —* pj decay rate at 
roughly the 5 — * 8% level, although the underlying effect has the potential 
to be larger. Contrary to other mechanisms of long distance physics which 
have been discussed in the literature, this yields a non-negligible modification 
of the B° —* p°7 channel and hence will provide an uncertainty in the ex- 
traction of Vtd- This mechanism also affects the isospin relation between the 
rates for B~ — > p~j and B° — > p°7 and may generate CP asymmetries at 
experimentally observable levels. 



I. INTRODUCTION 

The study of radiative rare decays of B mesons can provide valuable information on 
certain parameters in the CKM matrix. In particular, B — > p7 is thought to be an especially 
clean mode for extracting the Vtd matrix element. Although not yet observed, this mode 
should be accessible for study at B-factories. However, the smallness of the short- distance 
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B — > fry amplitude raises the concern that the experimental signal could be influenced by 
long-distance effects, with the possibility of their contributing a non-negligible fraction to 
the decay rate. 

The first studies @ of long distance contributions to rare B decays employed the vector 
meson dominance (VMD) model (cf Fig. 1(a)). As applied to the B — > /ry transition, @ 
this involves the decay B — > pV*° (V*° represents an off-shell neutral vector meson such 
as p°,u,(f), J/ip,- ■ •) with subsequent conversion of the vector meson to the photon. More 
recently, light-cone QCD sum-rules have been used to analyze the weak annihilation long- 
distance contribution to B — > fry. @ Here, the main effect arises from direct emission of 
the final state photon from the light spectator quark in the B meson followed by bq weak 
annihilation. 

Both these approaches suffer from some degree of theoretical uncertainty. Besides the 
usual model dependence of these methods, in both cases no allowance is made for final state 
interactions (FSI). However, the pj decay mode can be generated by the decay of a B meson 
into a hadronic final state which then rescatters into pj. It has recently been shown that 
soft-FSI effects are 0(1) in the large m;, limit and cannot legitimately be ignored. 

In this paper, we provide an estimate of a long distance component to B — > p7 which deals 
solely with on-shell transition amplitudes while explicitly taking the occurrence of FSI into 
account. Our mechanism can be viewed as a unitarity analysis, where a B decays into some 
intermediate state which then undergoes soft-rescattering into the final p7 configuration. Of 
course, we are not able to account for all possible intermediate states, so that in order to 
provide a rough estimate for the the effect we must analyse a few specific contributions. For 
definiteness, we consider the p°p° and p + p~ intermediate states in B° decay. In the language 
of Regge theory, fli the first of these proceeds by Pomeron exchange and is technically the 
dominant contribution, remaining nonzero in the heavy quark limit. However, the p + p~ 
contribution, whose rescattering is mediated by the p trajectory and is thus nonleading, 
can be numerically important because B° — > p + p~ is color-allowed (whereas B° — > p°p° 
is color-suppressed) and because the relatively low value of the B mass turns out to blur 
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somewhat the distinction between Regge-leading and nonleading contributions. Finally, 
since the corrections considered here are not themselves proportional to Vtd, their presence 
constitutes a potential source of serious error in the phenomenological extraction of Vtd- ® 
This underscores the importance of obtaining quantitative estimates of such effects. 



II. SOFT HADRONIC RESCATTERING 

Effects of final state interactions (e.g. Fig. 1(b)) are naturally described using the unitarity 
property of the 5-matrix, S^S = 1. This condition implies that the T- matrix, 5 = 1 + iT, 
obeys 
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ViscT B ^ f = - (f\T\B) - (f\T*\B) =^J2(f\Ti\I)(I\T\B). (1) 
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For the B — > p7 transition the contribution from the p intermediate state is the one which 
is most amenable to direct analysis, and we shall detail the B — > pp — > p7 component 
throughout this paper.Q This encompasses the p°p° and p + p~ intermediate states for B° 
decay and p + p° for B + decay. 

The final state interaction of Fig. 1(b) together with the unitarity condition implies a 
discontinuity relation for the invariant amplitude M.B^py, 

2Disc M B e ^ n e;(\ 3 y d (\ 4 ) = 

d 4 p 



E /^(2vrW-m 2 )5(( PB -p) 2 -m 2 ) 

x M B % p e;(A0e* (A 2 ) • e^e^Aa) M^tf <(A 3 )e;(A 4 ) , (2) 



Ai,A2 



x One should keep in mind, however, that for the S'-matrix to be unitary, inelastic effects such as 
diffractive dissociation (decay of the B into a p and a jet of particles which recombine into the 
final state) or any other intermediate state with suitable quantum numbers must also be present. @ 
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with rfW being the polarization tensors of p mesons in the loop. When expressed in terms 
of helicity amplitudes, the above discontinuity formula simplifies to 

IDisc A^ j b^ P7 (A 3 A 4 ) = 

£ / -%(27r) 2 5(p 2 - m 2 p )5{{ PB - P f - ml)M B -, pp (X 1 X 2 )M; p ^ (n (X 1 X 2 ; A 3 A 4 ) . (3) 
AiA 2 ^ > 

The FSI itself will arise from the scattering pp — > pry. We shall employ the Regge-pole 
description for this rescattering process, requiring both Pomeron and p-trajectory contribu- 
tions. 

To begin, however, we recall from Regge phenomenology the well-known invariant am- 
plitude for the scattering of particles with helicities {Aj}, 'i 

/ y-\ m/2 „-ina(t) if / „ \ a(t) 

where m = | A3 — A1I + IA4 — A2I and \J\ = 1. In the B — > pp weak decay, the pp state can exist 
in any of three helicity configurations, Ai = A2 = +1, Ai = A2 = —1 and Ai = A2 = (or in 

obvious notation ++, and 00). Moreover, the Pomeron (and near t — also the leading 

p) exchange does not change the helicities of the rescattering particles, i.e. Ai = A 3 , A 2 = A 4 . 
In view of this and noting that the photon helicity must have |A 7 | = 1, we omit the 00 helicity 
configuration from further consideration. This is equivalent to maintaining the condition 

of gauge invariance. Moreover, since parity invariance constrains the ++ and helicity 

contributions to be equal, we drop the A1A2 superscript hereafter and take for the residue 
couplings 7ll = 7^+ = 7. Throughout we use the linear trajectory forms, 

otp(t) = a° P + a'pt , ctp(t) = a° p + a' p t . (5) 

In the following discussion, we limit ourselves for clarity's sake to B° — > p°7 decay. The 
Pomeron part of the FSI occurs for the p°p° intermediate state. We obtain 

1 

Disc M^o^po^ = — — Mb^ p o p o / dt M* p o p o^ p o^ 
7p / s \ a ° p r ° 
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(—Y P M B o^ p a p0 [ dt e a 'p ln{s/so)t (cot[vra(t)/2] + i) , (6) 



where J = +1 for Pomeron exchange. Observe that the cotangent factor will diverge 
at some t = t such that a(t ) = 0. For the Pomeron, this occurs at t ~ —4.3 GeV~ 2 , 
which lies outside the forward diffraction peak. Such spurious behavior has been well-known 
since the early days of Regge phenomenology and is avoidable by restricting the range of 
integration to the diffraction peak, by employing a modified 'phenomenological' amplitude 
or by invoking daughter trajectories to cancel the divergence. We adopt the first of these 
procedures and find 

/ \ rv° -1 

,(P) 7p so ' ' x r ' 



Disc M%L P o, = H f l ^ (0 - 26 ~ °- 92 l) _ M B o^ p0p0 

1 1 327T a P m(s/so) \s J 



s x^p- 1 



= e P ( — ) M B o->po p o , (7) 

\Sq/ 

where s — 1 GeV. The numerical quantity ep is defined to encode the strength of the 
Pomeron-mediated rescattering. We discuss later how to determine the Pomeron residue 

function jp by fitting to experimental data. 

(p) 

The quantity Disc a4 b0 _^ must itself be inserted as input to a dispersion relation for 
o . Since the s-dependent factor in Eq. (0) is almost constant, the dispersion relation 
for the Pomeron contribution will require a single subtraction. Approximating a% ~ 1, we 
have 

= M^ pOj (0) + 6 fM B0 ^o p0 ln(l - £|) . (8) 

We now turn to the process B° — > p + p~ — > p°7. The formulas derived for the lead- 
ing Pomeron contribution of Eq. (0) are readily applicable to this case if one replaces the 
Pomeron trajectory by the p trajectory, but now with J = —\. This yields for the discon- 
tinuity function 

Disc Mf , = (-Y P M B o^ p+p - f° dt e a >^ (1 +itaa[7ra(f)/2]) , 

-(0.92-0.33 i) ( — ) M B o^ p0p0 



327r a' p \n(s/s ) \s 

( s \ a P _1 

e P — M B o^ p+p - , (9) 

\Sq/ 



where e p is analogous to the quantity ep of Eq. (f?|). The asymptotic behavior, 



Disc M { j';l_ tp0l /.s\""-' 



justifies in this case use of an unsubtracted dispersion relation for the amplitude, 



e p f°° ds e p y/s^ 2m p - m B 

7T JAm 2 p \fS{S — m B ) 7T VTLp 2m p + VTLb 



(10) 



where we approximate a° p ~ 0.5. The above discussion is extendable in like manner to 
B + — > p + 7 decay. 



III. THE WEAK DECAY VERTEX 



The most general form for the weak amplitude B — > pp is 

T h 

^B(p)^ p ( fcl )p(fe 2 ) = e tl {k 2 )e* v {k 1 



b c 
ag P u H TyV»Pv + i—zCiivapkfpP 



(12) 



m P m p 

The quantities a, b, c can be interpreted as partial-wave amplitudes as they exhibit the 
respective threshold behavior of S, D, P- waves. To proceed further requires knowledge of 
{a, b, c} for both the B° — > p°p° and B° — > p + p~ weak decays. Since no data yet exists for 
the B — > pp transitions, we must determine {a, b, c} theoretically. We have employed the 
BSW description of the nonleptonic B transitions. ® In this model, the amplitude for the 
B°p + p~ transition is 
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( ' F VubV* d ai (p-(k 2 )\dYu\0)(p + (k 1 )\u lp (l + l5 )b\B°(p)) , (13) 



lBO(p)-*p+(fe 1 )p-(fc 2 ) - 

where a\ ~ 1.03. @ We adopt standard notation for the matrix element 



(p-(k 2 )\uYd\0} = f p m p e^(k 2 

2i 



(p+ih^blB^p)} 



tub + m. 



(14) 



(p+ih^^blB^p)) = (m B + m p yjk 1 )A 1 (k 



6*{h 



2e*(fcQ ■ k 2 

Jv 1 ) 



k 2p m p 
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We + rn p 
A 3 (k 2 2 ) - A (k 2 2 ) 



In what follows we assume the k\ behavior of the form factors to be of the simple pole form, 



and for k\ = m 2 p we obtain 



V(k 



2\ hv 



2! 1 - kl/ml, 



Mkl) = , ^ 2 -h Ai (< = 0,..., 3), (15) 
with hy ~ 0.33, h Al ~ /iy^ ~ 0.28. § This gives for the amplitudes in Eq. ([12|) , 

a ~ -/|K6K*d m P a i/p( m s + m p )h Al , 

& - -%KbK>i/p — 2 — £ — h A 2 , (16) 



V2 U ° ud lJp m B + m 



c ^ -^VubV^aJp — fry . 

The set of relevant helicity amplitudes can be written using explicit form of the polarization 
vectors. The connection between the transverse helicity amplitudes and form factors in 



Eq. (|TJ) is 



M++ = a + Vx^l c , M~; = a- Vx^l c (17) 



where x = {m%— 2m 2 ) /2m 2 . For completeness, we note that branching ratios for the B —>■ pp 
decays calculated from Eq. ([17]) amount to B^o^ p o p o = 4.8 • 10~ 7 , B§o^ p + p - = 2.1 • 10~ 5 and 
B^-^p-pO = 1.5 • 10~ 5 , where all three pp helicity configurations have been summed over. 

IV. NUMERICAL RESULTS 

Before numerically estimating the long-distance effects in B —>■ p~f, we first recall the de- 
termination of the short-distance contribution which arises from the effective Hamiltonian B 

G 



n s - d - = -^(v tb v; d cf o 7 + 



7 =^dap U (l- l5 )F^b , (18) 



rr 

where Cj = —0.306. The amplitude for the short-distance B — > /ry contribution is 



^V tb V t * d C™ ^| e*»(q, X^ih, X p ) 2F b (0) [e^k^ - % (m B E,g, u - p p q v )\ , (19) 

where A 7 , X p are the respective photon, rho heli cities and F s (q 2 = 0) is a form factor related 
to the B-to-p matrix element of O-j. It is estimated from QCD sum rules @ that F s (0) ~ 0.17 
for B° — > p°7 and _F S (0) ~ 0.17-\/2 for B~ — ► p~7- This i? — > p7 matrix element can be 
written in a form similar to Eq. (|l^) with the identification 

« s ' d - = i^V tb V t * d cf e £±m B E,F s (0) 

c sA - = ~b sA - = i^my th V; d Cf e £^F s (0) . (20) 

These can be related to the set of helicity amplitudes by using Eq. ([17]) which then contribute 
to the decay rate for B — > /ry as 

rw- d - = ^|(l^ + l 2 + l^:l 2 ) , (2i) 

where q is the photon momentum. 

In our numerical work, we have adopted for the CKM matrix elements the values 



\V cb \ = 0.0393 , IVJ = 0.08 \V cb \ , 1.4 < 



td 



< 4.6 . (22) 



V ub 

The biggest source of uncertainty for the short- distance contribution is clearly the magnitude 
of Vtd- The above range leads to a range of branching ratios for the short distance component 

^V 7 = T -%^ - 10- 7 — 10- 6 , (23) 

1 BO 

with i3^ d ^ p _ 7 = 2£>£ol^ p o 7 - It is this large theoretical spread which has motivated an 
experimental determination of £>b^ P7 as perhaps the best solution to the Vtd problem. At 
present, however, even the largest theoretical values are still considerably smaller than the 
experimental upper bounds Sbo_ >(0 o 7 < 3.9 x 10~ 5 and Bb-^ p - 1 < 1-1 x 10~ 5 . @ 
Upon including both short-distance and long-distance contributions, we obtain 

Vb ^ = ^(\ M tl + M + d + l 2 + \m- a . + mxI 2 ) . (24) 



The long-distance amplitudes are in turn summed over contributions from the Pomeron and 
p trajectories, 

Mt d X2 =P {P) M x p $ +P [p) MX 2 - ■ (25) 

Using Eqs. (8), (11) and dropping the subtraction constant but keeping the terms which are 
clearly related to rescattering, we obtain 

0(P) = ££ Ul - ^) , 0« = ^ In 2m " ~ mB . (26) 
tt v Amy 7r m,B 2m p + mg 

Also needed for estimation of the FSI effect are the residue functions jp and j p . For the 
Pomeron case, we use experimental data on photoproduction of p° mesons from a nucleon 
target. Upon calculating the differential cross section da/dt at t ~ from the matrix 
element of Eq. (Q) and applying the quark- counting rule to relate Pomeron-pp and Pomeron- 
pp couplings, we obtain 7p ~ 4.53. Actually we use twice this value as the final-state photon 
can arise from either scattering vertex. For p exchange, we use 7p — > p + n data together 
with the isospin relation 7 P (7p° — > P + /l P {.lV — ► p + n) ~ 2 to obtain 7 P ~ 7.6. 



The recent review of the CKM matrix given in Ref. JTT[] cites p = 0.05 and rj = 0.36 as 



providing the best fit to current data, where p and r\ are the Wolfenstein parameters 

V ub = AX 3 (p - irj) , V td = A\\l-p-iri) , (27) 

and which corresponds to |Vt<f/Ki&| — 2.8. Assuming these central values, we obtain for the 
short distance and short-plus-long distance branching ratios the respective values 

BSoUpo 7 - 7.3 x 10~ 7 , B^ p0j ~ 7.8 x 10" 7 . (28) 

Figure 2 displays the effect of varying the p, rj values over their allowed range . The upper 
and lower curves define the band of allowed values for the ratio {13^ ~— B%f ) / i3^oi^ p o 7 - 
The upper curve is seen to be about 8% for virtually the entire physical range of p. 
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V. CONCLUDING REMARKS 



Final state rescattering effects may sometimes, although not always, modify the usual anal- 
ysis of decay processes. The situations where rescattering can be important are those where 
there is a copiously produced final state which can rescatter to produce the decay mode 
being studied. In our case the most important channel is B — * p + p~ , which is color allowed. 
It is also proportional to V u b instead of Vtd so that it is clearly a distinct contribution, and 
one that becomes more important if Vtd is at the lower end of its allowed range. Our results 
should also be adjusted upward or downward if the measured branching ratio proves to be 
larger or smaller than that predicted above using the BSW model. 

The magnitude of the pp — > p7 soft-scattering for B° — > p°7, summarized in Eq. (f28[), 
is seen to occur at about the 8% level. This differs from the QCD sum rule estimates of 
Ref. || which found a 10% long distance contribution to Bb+-> p +~, but only a 1% effect for 
Bsa^pOry. Our result is especially noteworthy in view of the many other possible contributions 
to the unitarity sum and the finding of Ref. || that multi-particle intermediate states are 
likely to dominate. Thus, although it is not possible at this time for anyone to completely 
analyze the long distance component, it seems plausible that the FSI contribution has the 
potential to occur at the 10% level and perhaps even higher .0 We concur with comments in 
the literature H' @> II that a deviation from isospin relations based on the short distance 
amplitude {e.g. B B -^ p -^ = 223 B o_ j . p o 7 ) will be evidence for a long distance component, and 
such is the case here. Defining a parameter A to measure the isospin violation as 

A^l-i^O: , (29) 

we display in Fig. 3 the dependence of the effect upon the Wolfenstein parameter p. The 
isospin violation is found to be largest for large p. 

2 However, we do not claim the dominance of the long-distance contribution over the short-distance 
amplitude, as was done in Ref. @ for B — ► 7T7T. 
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The FSI mechanism described here will not markedly affect the B — > K*^ rate. The 
leading CKM contribution would involve a D* — D* intermediate state and thus be very 
suppressed since the soft-scattering requires a Reggeon carrying the quantum numbers of 
the D meson. The non-leading CKM amplitudes, involving contributions from light-quark 
intermediate states, are proportional to V u b and are likewise very suppressed. 

Finally, we point out the implications of our work to obtaining a CP-violating signal. 
Since B — > pp decay is governed by CKM matrix elements which differ from those describing 
the short distance B — > p7 transition, the necessary condition for CP violation is satisfied. 
We consider the CP asymmetry 

a = l r B->P7 ~ r B-^o 7 l = \Bb-> p ~/ - ffg^ P7 | ^ 

We have studied the resulting effect numerically. The FSI are found to increase the amplitude 
for B° — > p°7 but reduce it for B° — > p 0r y. This is what gives rise to the asymmetry and 
we find a ~ 7%. As with our other results, we take this magnitude as an indication that 
experimentally interesting signals might well exist and should not be ignored in planning for 
future studies. 

We would like to thank Joao M. Soares for reading the manuscript and useful comments. 
This work was supported in part by the U.S. National Science Foundation. 
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FIGURES 
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Fig. 1. (a) VMD mechanism and (b) FSI mechanism. 
FIG. 1. B — ► p7 via (a) Vector dominance process and (b) final state interaction effect. 
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FIG. 2. Dependence of the FSI effect upon the CKM parameter p. 
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Isospin Violation, % 




FIG. 3. Dependence of isospin violation upon the CKM parameter p. 
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